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Abstract 



The magnetic inverse source problem of reconstructing the positions and cur- 
rents of very long parallel conductors is considered in a two-dimensional situation, 
with applications to power line measurements. The input data is the magnetic field 
on a contour surrounding the conductors to be reconstructed. Using a scalar-vector 
Green identity, an explicit reconstruction algorithm is derived. The numerical im- 
plementation of the algorithm is described and simulation results are presented, 
demonstrating the influences from numerical errors and uncertainties in measure- 
ment data. The algorithm can handle an arbitrary number of conductors, but sta- 
bility problems associated with the illposedness accelerate with increasing number 
of conductors. Mathematically, the Green identity approach removes the influence 
of external disturbances and thus have potential usefulness in current reconstruction 
for determining optimal sensor positions and how to process measurement data. 

Keywords: current reconstruction; harmonic function; magnetic inverse source prob- 
lem; power line measurements; scalar-vector Green identity 

1 Introduction 

In electric power engineering, magnetic inverse source problems appear in contact-free 
measurements of currents in power lines [Il[2], bus bars [3H5] and cables [BHSj. The 
magnetic field is measured by a set of sensors located in the vicinity of the conductors, 
and the currents are determined by inversion of the direct map from the currents to the 
magnetic field, using a suitable model of the problem. Due to the illposedness of such 
inverse problems, the results become highly sensitive to errors in the measured magnetic 
field and to deficiencies in the models. Hence, much effort has been spent on reducing the 
impact of disturbances from exterior magnetic fields [3l[5] , and on finding optimal sensor 
positions [HITj. In electric power engineering, the focus has been on reconstructing the 
currents in individual conductors. Owing to the fact that in many situations the conductor 
positions are known or determinable by other means than magnetic field measurements, 
the simultaneous reconstruction of positions and currents has devoted less interest. An 
exception is the early work in [2], where also the conductor positions are reconstructed, 
from experimental data. 

Another application of magnetic inverse source problems is in biomagnetics. For example, 
in magnetoencephalography (MEG) epileptic activity is detected as current sources within 
the human brain [9]. With MEG being a tool providing input data to radiotherapy or 
surgical treatment, usually the localization of the current source is more important than 
the detailed determination of the current. In e.g. MEG, the localization problem can 
sometimes be solved explicitly by means of methods based on using harmonic functions 
in combination with Green identities [T014T2] . Such explicit methods typically require 
measurement data on closed surfaces that completely surround the region containing the 
sources. In MEG, magnetic field data can be obtained over a quite large solid angle, which 
together with data-continuation facilitates explicit identification methods [12]. However, 
an obstacle is that the current sources are submerged into a conducting medium while the 
magnetic field is measured in another, non-conducting, region that typically is air. Hence, 
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for quasistatic fields, measurement data cannot be continued into tfie source region, in a 
stable manner. One way around that problem is to use a simplistic model (e.g. model the 
brain as a sphere) that allows the configuration to be replaced by an equivalent source in 
free space [13]. However, such an approach may introduce large modelling errors. 

In many electric power applications, the conductors are located in free space, i.e. the same 
region in which the surrounding magnetic field is measured. Thus, one will not encounter 
the difficulties with data continuation that appears in MEG problems. Hence, it is of 
interest to investigate whether the explicit localization methods, considered previously in 
MEG, can be used in power engineering for reconstructing currents and localizations of 
conductors, in overhead power lines, bus bars etc. 

In this paper, we consider the two-dimensional approximation of the magnetic inverse 
source problem of reconstructing the positions and currents of very long parallel conductor 
structures. In section O we formulate the problem and present the theory behind the 
explicit method. The reconstruction algorithm is described in Section [3l The numerical 
implementation and simulations results are presented in Section HI and Section [5] contains 
the conclusions. 



2 Problem formulation and preliminaries 

We consider a two-dimensional approximation of a parallel conductor structure. The cur- 
rents distributed over the cross section of the structure are modelled as a two-dimensional 
longitudinal current density 

J{r) = J{x,y)z (1) 

where r = xx + yy. x, y and z denote the unit vectors in the cartesian coordinates; 
see Figure [1] The current density J may consist of several disjoint parts and the inverse 
problem is to reconstruct the part internal to a measurement contour C by using magnetic 
field data on C. 




(a) General longitudinal current density (b) Line-current model of the current 

in the reconstruction region. density in the reconstruction region. 



Figure 1: The problem geometry. The inverse problem is to reconstruct the current 
density inside the contour C using magnetic field data on C. S is the surface bounded by 
C, fi and t are the unit normal and unit tangent, respectively, on C. 



3 



In the reconstruction region S, the surface bounded by C, we have no access to magnetic 
field data. It is assumed that C passes between disjoint parts of J, which means that on 
and in the neighborhood of C we have J = 0. In practice, the current density exterior to 
C represent nearby parallel conductors and/or the infiuence from a parallel ground. 

The preliminary analysis is carried out assuming that J has a general distribution in 



the reconstruction region; see Figure 1(a) The reconstruction algorithm, described in 



Section [3], is restricted to the case when J is a finite number of line-currents; see Figure 



1(b) Assuming slowly varying currents, we neglect propagation effects and apply static 



analysis methods. Hence, the magnetic field B fulfils 

V ■ B{r) = (2) 
V X B{r) = /io J(r) (3) 

where V = xd/dx + yd/dy. From ([T]), and it follows that J generates a transversal 
magnetic field: B = B^x + Byy. If we measure amplitudes and phases of time-harmonic 
fields, instead of instantaneous values, J and B represent complex phasors, with the 
time- dependence exp(jcjt) assumed and suppressed. 

2.1 Integral equation for J 

Consider the following Green identity for a scalar and a vector function [T3] : 

j [GV^B - BV^G] ds 

= ^ [G {V ■ B) h - {h ■ B)VG + G {V X B) X h - (h X B) X VG] dl (4) 

where the scalar function G is a function to be specified and fi is the from S outwardly 
directed unit normal on C; see Figure [H (jl]) holds generally, and with no current density 
on C ([3]) implies V x S = on C Thus, (jl]) reduces to 



J [GW^B - BV'^G] ds = - jf [(n ■ B) VG + (n x B) x VG] dl 



(5) 



which means that on C we do not need to measure derivatives of the magnetic field. 
However, over S the magnetic field is not accessible for measurements, and to handle that 
in ([5]) we require the function G to be harmonic over S, i.e. 

V^G = (6) 

Hence, ([5]) simplifies to 

/ GV^Bds = - (f [{n ■ B)VG + {fi X B) X VG] dl 
Js Jc 

From ([2]) and ([3]) it follows that the Poisson equation for B becomes 
V^S = V (V ■ B) - V X (V X B) = -/ioV X J 



(7) 
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which, by integration by parts over S utihzing that J = on C, yields that the left hand 
side of 

GV^Bds = -/io / GV X Jds = -i_iQ [ Jx VGds (9) 
Js Js 

With J xVG = Jz xVG and -z x{z x VG) = VG, we obtain using ([7]) and that 

JVGds = -—z X (f[{h- B) VG+{hx B)x VG] d/ (10) 
/^o Jc 

Considering the right hand side as known, ffTOl) is an integral equation for the current 
density J(r), where the gradient of the harmonic function G(r) is the kernel. 



2.2 The choice of the harmonic function G 

Aiming at an explicit reconstruction algorithm, the harmonic function G(r) in ffTOl) must 
be reasonably simple. Similar to in [IHIIIB], G is defined as 

G{r) = G{w{r)) , (11) 
w = k-r (12) 

where the complex vector k = k^x + kyi; fulfils k ■ k = 0. Denoting dG/dw = f{w), we 
obtain 

VG(r) = f{w) Vw{r) = f{w) k (13) 

V^G(r) = V ■ VG(r) = k ■ V/(u;(r)) = ^k ■ Vw{r) = ^k ■k = (14) 

Hence, for ([6]) to hold, the vector k must be complex. For three dimensions, the general 
form of k is given in [10]. For two dimensions, we conclude that the general form of k is 

k = a{x + isy) (15) 

where a G C is an arbitrary constant and s = ±1. Using f[T^ in fllU]) . it follows that 

k I Jfds = -—z X (p [{h- B)k + {h X B) X k] fdl (16) 
Js /^o Jc 

Since the two components of (|T6l) are in effect two identical equations, we dot-multiply 
( !T6|) with k* {* denotes the complex conjugate) which after some algebra yields 

! Jfds = — (f [isn- B + i- B] fdl (17) 
Js /^o Jc 

Since we are working in two dimensions, it follows from ( IT2|) and ( IT5|) that in (ITT]) 

w = a {x + jsy) (18) 

which is an analytic function of either of the complex variables x + iy or x — iy. Hence, 
([6]) becomes fulfilled if G{w) is an analytic function for all w corresponding to r G 5. 
Consequently, f{w) = dG/dw is also an analytic function. 

Note that with / = constant, the property n ■ Bdl = V ■ Bds = yields that f|T71) 
reduces to the well-known circulation law B ■ tdl = /io IgJ' ^ds. Hence, in this two- 
dimensional situation, f lT7|l can be considered as a generalized circulation law, through 
which more information about the enclosed current density can be obtained by altering 
the harmonic function /. 
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3 Algorithm for reconstructing line-currents 



Due to the existence of so-called silent or non-radiating sources, which produce zero fields 
outside their region of support, generally formulated inverse source problems typically ex- 
hibit non-unique solutions; a property utilized practically in so-called magnetic signature 
reduction [15]. In Appendix El we give examples of silent current densities that cannot be 
determined using the integral equation (fTOj) . One way to achieve uniqueness is to shrink 
the space of possible solutions. 

From now on, we restrict the interior current density J{r) to be a finite number of 
line- currents: 



N 



J{r) = Y,InKr-rn) 



(19) 



n=l 



where {In}n=i currents and {rn}n=i conductor positions; see Figure 1(b) 

(6() is the Dirac-delta function). The inverse problem thus becomes to reconstruct the 
currents and the conductor positions. In this study, we assume that is known; for the 
determination of A^, see [TT]. 



(20) 
(21) 



With J = Jz, we obtain from (jH]) and ( fT9i) that in the reconstruction region S 



N 
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V^5^ = -/io ^/„— 6(r - r„) 

n=l ^ 



N 



V25, = /io5^/„g^6(r-r 



n=l 



Hence, mathematically, the interior source of the field components is a finite number of 
dipolar sources, and in [11] it has been shown that this kind of inverse problem has a 
unique solution. 

The reconstruction algorithm has originally been derived in [11]. Here, we give an alter- 
native derivation, adapted to the present context. With f{w) being an analytic function 
it follows that [/(w)]™ , m = 0, 1, 2, . . . are analytic functions as well. Let = [/(w^)]™, 
where Wn = k ■ r„. Hence, we obtain from f[T7|) and f[T^ that 



N 

Y,Un = -f[jsh.B + i.B] rdi 



m 



0,1,2, 



(22) 



where bm are constants determined from the magnetic field measurement on C and the 
choice of f{w{r)). In (l22l) . we have 2A^ unknowns in {In, fn}n=i- {fn}n=i have been 
determined, the currents can be obtained from a linear equation system generated by A^ 
different values on m. For example, using m = M, M + 1,...,M + N — 1 (where M > 0), 
we introduce the following matrix and vectors: 







fM ' 

J N 




'h 


Fm — 






1 = 






rN-l+M 

Jl 


rN-l+M 

■ Jn 




In 



'M 



^N-l+M 



(23) 
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Hence, if F m is known the current vector / is determined from the equation 

FmI = bM 



The complex numbers {fn}n=i are the roots of a polynomial equation: 



Introducing the vector c containing the polynomial coefficients: 



(24) 



(25) 



Co 



CN-1 



(26) 



we obtain using ( 125|) that 



c^Fm 



(T denotes matrix transpose). Hence, using (124|) and (p7|) . the scalar 

r rN+M rN+An J _ 1 



(27) 



(28) 



Now, using different values of M in the relation bj^c = —Bm+n, we obtain a linear 
system from which the coefficient vector c can be determined. For example, using M = 
L, . . . , L + N — 1 (where L > 0), we construct the matrix 



"L+N-1 



(29) 



in which way the coefficient vector c follows from the equation 

Clc = -bi+N 

With c determined, {fn}n=i follows from fl2^ . whereafter {In}n=i follows from 



(30) 



4 Reconstruction results 



4.1 Numerical approximation of the measurement integral and 
reconstructions from clean data 

The input to the reconstruction algorithm is the 6m-coefficients, obtained by integration 
over the measurement contour C (see f lT7|) ): 



bm = — (f[isn-B + i-B] rdl, 
f^o Jc 



m = 0,1,2,... 



(31) 
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(a) Measurement contour and con- (b) A straight segment AC of C be- 

ductor locations used in the numer- tween two measurement points, 

ical examples. 



Figure 2 



In a practical situation, magnetic field data is obtained only at a finite number of points, 
wherefore (13T|) must be evaluated approximately. Here, we choose to distribute the mea- 



surement points uniformly around a circle with the radius Rmcas', see Figure 2{a] 



Since the integral in fl3T|) is independent of the curve C as long as it circumscribes only the 
currents to be reconstructed, we chose C as straight segments between the measurement 
points. On a segment AC between points ri and r2, see Figure 2(b), we obtain 

= ri (1 - t) + rat, d/ = |r2-ri|dt (32) 

where the parameter t G [0, 1]. At ri and r2, the measured magnetic fields are denoted 
Bi and B2, respectively. On AC the magnetic field is interpolated linearly: 

B{t) = Bi{l-t) + Bat, r(t) G AC (33) 

Utilizing that h = i x z, the contribution to (I3T]) from AC becomes 

A6„ = / [B,{t) (t,. + jst,) + By{t) {ty-isQ] [f{r{t))rdt (34) 

fJ'O Jo 

Even though many choices of /(t) admit explicit evaluation of the integral in f lM|) . we 
find it from the versatility point of view better to use a numerical quadrature routine; the 
numerical accuracy is anyhow superior to the approximation in (l33l) . 

The analytic function f{w) = f{x + jy), we chose as 

f{x,y) = expl^-^ — -] (35) 

Hence, once the evaluations /„ at the conductor positions have been found, the conductor 
coordinates are obtained as 

Xn = ^measlm {lu /„} (36) 
Hn = --RmeasRe {lu /„} (37) 

Note that, in addition to being analytic, the function f{w) must be invertible in the 
reconstruction region S. 
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As the test case, we consider three internal conductors, whose positions and currents 
are to be reconstructed, and two external conductors that act as disturbance sources. 
The coordinates and phasor currents of the conductors are given in Table [1] (conductor 



positions also illustrated in Figure 2(a)). 



Table 1: Locations and currents for the conductors used in the numerical examples. 



Conductor data 


Inside C 


Outside C 


Coordinate x/i?mcas 


-0.5 





0.5 


-1.5 


1 


Coordinate y/Rmcas 


-0.5 


-0.5 


-0.5 


-0.5 


1 


Current I/A 


-j 


2 


1 


-1 


2j 



For the generation of artificial magnetic field data, we use the formula 

rn) 



N+Ncyct T - ( 



271 



n=l 



r„ 



(38) 



To check the accuracy of the linear interpolation, we take the input data from Table [T] 
and use fl38|l in fl33|) and flMl) . and compare with the exact results, obtained by using the 
input data in the left hand side of (1221) . For a smooth magnetic field, the dominant error 
in the linear interpolation f l5^ is quadratic in the segment length, wherefore the rate of 
convergence of the integral for 6^ is expected to be proportional to the inverse square of 
the number of measurement points; this was also supported by observations. Hence, we 
use Richardson extrapolation to improve the numerical accuracy. We sample the circular 
measurement curve at every 10°, as depicted in Figure 2(a) , and let b"^ denote the results 
when using all 36 measurement points. Then we use data from 18 points at even and odd 
multiples of 10° and denote the corresponding results fe^*"^ and b'^'^, respectively. The 
average of the extrapolations from the even and odd sets becomes 



lextrapol 



b^! - bl 



Lodd 



6 



(39) 



In Table |2l we compare the results, when using m G {1,2,3,4,5,6}; the necessary and 
sufficient number of m- values for reconstructing three line-currents. We conclude that the 
Richardson extrapolation improves significantly the numerical accuracy. 

To evaluate the stability of the reconstruction algorithm against the numerical approxi- 
mations, we first test the algorithm using clean measurement data. Referring to Section 
[31 we use the values A^ = 3,M = L = 1. We chose three different values of the number, 
Nraca.si of measurement points. For the internal conductors, the results for the positions 
are presented in Table [3] and the results for the currents are presented in Table HI We 
observe that each doubling of the number of measurement points reduces effectively the 
errors in the reconstructions. 



4.2 Reconstructions from noisy data 

To evaluate the stability of the reconstruction algorithm against random measurement 
errors, we add artificial noise to the clean data generated by (1381) . On the real and 
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Table 2: Results, in the unit A, of extrapolations of the coefficients bm, evaluated by flM|) . 
b"^^^ are from measurement points at {0°, 20°, . . . 340°}, on the circular measurement curve 
C. b^*^ are from measurement points at {10°, 30°, . . . 350°}. b'^ are from all measurement 
points, fj^trapoi ^Yie extrapolations using f l5U]l . b"^^^^ are the exact results, calculated 
directly from the conductor currents, m is the power of the function /. 



m 


ijCvcn 

■m 


ijOdd 

■m 


^all 

■m 


ijCxtrapol 




1 


1.027- 2.132j 


1.027- 2. 133j 


1.052-2.211j 


1.060 -2.237j 


1.060- 2. 237j 


2 


1.494- 3.575j 


1.494 -3.573j 


1.630- 3.711j 


1.675- 3.757j 


1.681- 3. 756j 


3 


3.302 -4.793j 


3.312 -4.779j 


3.937- 4.797j 


4.147- 4.801j 


4.176- 4. 787j 


4 


8.215 -4.648j 


8.248 -4.612j 


10.331 -3.934j 


11.031- 3. 702j 


11.134- 3.644j 


5 


19.169- 1.171j 


19.238- 1.112j 


24.727+1. 45 7j 


26.568 + 2.324j 


26.834 + 2.469j 


6 


40.514+8.492j 


40.608+8.565j 


52.670+14.760j 


56.706+16.837j 


57.221+17.050j 



Table 3: Reconstruction of conductor positions when using clean magnetic field data, 
from 72, 36 and 18 measurement points, respectively. The results are presented as dis- 
placements from the true values, in units of i?m, the radius of the measurement curve. 





•^true/ Rmeas 


2/truo/ Rmcas 


-0.5 0.5 


-0.5 -0.5 -0.5 


' meas 


Ax in % of Rmcas 


Ay in % of i^^cas 


72 


0.01 


0.04 


-0.02 


0.05 


-0.03 


0.04 


36 


0.16 


0.60 


-0.24 


0.65 


-0.40 


0.64 


18 


1.79 


4.25 


-1.30 


6.08 


-3.52 


5.91 



imaginary parts of both magnetic field components, B^j. and By, we add Gaussian noise 
with zero mean value and with standard deviation determined as 

Cr = arcf( I -Bmcas I) (40) 

In fHU]) . (Jrcf is a reference parameter and {\B^cas\) is the average over the measurement 
points of the magnitude of the magnetic field. 

In the reconstruction examples, we consider two different levels of noise: (Xref = 0.01 and 
(Trcf = 0.05, and for each level we test with A'mcas = 18 and A'meas = 72 measurement 
points, respectively. To investigate the spread in the results we have for each case made 
50 simulations. The results for the conductor positions are presented in Figure [3] and 
the results for the complex phasor currents are presented in Figure |H We see that in- 
creased noise levels increase the spread in the reconstructions. Increasing the number 
of measurement points has less impact on reducing the spread, but affects strongly the 
average positions, which tend to coincide with the results obtained when reconstructing 
from clean data; see Table [31 

5 Discussion 

Although the restriction to a line-current model has removed the non-uniqueness, the 
inverse problem is still illposed, in the sense that the algorithm is sensitive to both nu- 
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Table 4: Reconstruction results for conductor currents when using clean magnetic field 
data, from 72, 36 and 18 measurement points, respectively. 



Currents/ A 


/i/A 


/2/A 


/3/A 


True values 


0-j 


2+Oj 


-1+Oj 


TV — 72 

^ ' meas ' ^ 


0.0019 - 1.0008j 


1.9993 - 0.0009j 


-1.0012 + 0.0017j 


-Ameas — 36 


0.0274 - l.OOQlj 


1.9874 - 0.0149j 


-1.0145 + 0.0243j 


A^meas = 18 


0.2489 - 1.0419j 


1.7706 - 0.1944j 


-1.0187 + 0.2384j 



y/R. 




-0.7 



Figure 3: Reconstructions of the positions of the internal conductors using data con- 
taminated with Gaussian noise. In each case we show 50 reconstructions. Green boxes: 
A^mcas = 18 and a^ci = 0.05. Blue boxes A^meas = 18 and a^ei = 0.01. Black diamonds 
A'^mcas = 72 and a^d = 0.05. Red diamonds A^meas = 72 and cxref = 0.01. 



merical errors and uncertainties in the measured data. Increasing the number of conduc- 
tors, beyond three conductors considered in the numerical examples, the stability problem 
associated with the illposedness increases with an accelerating rate. Hence, if positions 
and currents are to be reconstructed simultaneously, a very large number of measurement 
points will be needed and the magnetic field sensors must have a very high precision; 
similar conclusions are drawn in |2]. 

The reconstruction algorithm has certain similarities with the algorithm used in [2], al- 
though their algorithm is based directly on point measurements. In our case, the con- 
struction of the algorithm from a Green identity yields that the integral in the right hand 
side of f|T7|) has the property that the part of the magnetic field originating from currents 
outside the reconstruction domain gives no contribution; mathematically, our algorithm 
complete cancels the influence of disturbance sources. 

As mentioned, in most electric power applications, the conductor positions are known 
or determinable by other means. Hence, despite difficulties with the full reconstruction 
problem, our algorithm can be applied to the reduced problem of reconstructing the 
currents only. In that case, the coefficients in f l22|) are evaluated from the known 
conductor positions, and the currents are determined directly from (^^. Thus, since 
the mathematical formulation cancels external disturbances, the measurement integral in 
f l22]) can be used as guideline for determining optimal sensors positions [HEIIT] and for the 
processing of measurement data. 
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1 I r 




Re{J}/A 

Figure 4: Reconstructions of tlie pliasor currents in the internal conductors using data 
contaminated with Gaussian noise. In each case we show 50 reconstructions. Green boxes: 
-^mcas = 18 and (Trcf = 0.05. Blue boxes A^'^eas = 18 and a^-cf = 0.01. Black diamonds 
Nmeas = 72 and (Trcf = 0.05. Red diamonds iVmeas = 72 and a^d = 0.01. 
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A On the non-uniqueness of the current density 

Here we give two examples of when the current density J cannot be reconstructed uniquely 
using (fTOj) . 

A.l Surface currents 

Let So be a subsurface, completely interior to S, bounded by the contour Cq. We will 
show that on Co one can always put a surface current density Ko = KoZ that exterior 
to Co completely cancels the magnetic field from the current density Jo = JqZ in So- 
From ([2]), the transversal magnetic field Bq has the vector potential Aq = AqZ. Hence, 
Bo = V X Aq = VAq X z, and using ([3]) it follows that 

VMo = -fXoJo, re So (41) 

In consistence with no external field and absence of magnetic surface charge densities, we 
have on Cq that ho ■ Bq = ho ■ (VAq x z) = {z x ho) ■ VAq = to ■ VAq = 0. Hence, 

Ao = do = constant, r E Co (42) 
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(HTj) and fl42p are recognized as the interior Dirichlet problem for Poisson's equation. 
Hence, Aq is uniquely determined in Sq expect for the undeterminable additive constant 
ao in fj^2|) . Since Oq does not affect the magnetic field, we can set ao = 0. With the 
magnetic field determined uniquely, the required surface current density on Co, to achieve 
zero external field, follows from the tangential magnetic boundary condition: 

1 11 z 
Ko = —no X (0 - Bo) = no x Bo = no x (VAo x i;) = — (no ■ VAq) (43) 

or 

Ko = —no ■ VAo (44) 

Hence, for an arbitrary choice of current density Jo in '^o (which may also include surface 
and line currents) we always have a unique surface current density Kq on Co such that Jo 
and Ko produce no external field. This result is a two-dimensional analog of the magnetic 
signature reduction problem that has been discussed in [T3] . 

Since Jo and Kq produce a zero field on C, i.e. in the right hand side of fllUj) . their 
contribution in the left hand side of ffTOj) must vanish, i.e. 

f Jo VG ds+ ([ Jio VG dl = (45) 

-So J Co 

Expressed in index notation, (14T]) and (144|) . together with the Gauss theorem, imply that 
becomes proportional to 



dG d dAo , f dGdA, 



5q dxi dxj dxj Jc,, dxi dxj 



ds + (b -r—— — -Ui dl 



1 5^ dxi dxj dxj J 5p dxj \ dxi dxj J J^^ dxj dxj dxi 
Integration by parts, with Ao = on Co and using (|6]), yields 

f dAo d dG ^ f , d d dG ^ f , d , 

Js^^ dxj OXj dxi oxj OXj OXi oxi 

which verifies (HSl). 



A. 2 Volume currents 

If we have a coaxial current density Jo(p), with support within the distance ao from a 
symmetry axis at p = 0, that transports zero total current, i.e. fulfils the condition 

ao 

Jo(p)pdp = (48) 

the magnetic field Bq = when p > ao- Hence, a superposition of such current densities, 
having parallel symmetry axes, produces no magnetic field outside its overall region of 
support. As a consistency check, we verify that a silent current density subject to 
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yields a zero contribution in the left hand side of fllOp . Obeying the Laplace equation, G 
has in the region p < the following general solution, in the polar coordinates {p, 

G(p,^) = ^^?„pHej-^ (49) 

meZ 

where Qm are coefficients determined by the specific choice of G. Hence, 

VG = + = V (|m| p + ]m(p) pH-iej-v' (50) 

Consequently, we must verify that 

JoVGds = Y,9m Jo(p)p''"I^Vdp / (|m|p + jm<^)e5'"^d<^ = (51) 
Js ^^-T, Jo Jo 



Expressed in the fixed cartesian basis {a;, y}, the unit vectors {p, c^} are functions of e^^'''. 
Hence, for m 7^ ±1 the angular integral vanishes, while for m = ±1 the condition (148!) 
applies to the radial integral, which verifies (15 ip . 
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